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Abstract— Indonesia is a country with the highest malaria cases in the world. North Sulawesi is known as 
one of the malaria endemic areas in Indonesia. Malaria can only be transmitted through the bite of the 
Anopheles sp. Thus, the high case of malaria infection in an area is linear with the high population of 
Anopheles sp. The identification method to the species level that has high accuracy is by molecular 
identification using the cytochrome oxidase sub unit 1 (COI) gene. Based on the COI gene, the 
mitochondrial DNA of Anopheles sp from Tombatu was 92% similar to Anopheles maculatus 
[KT382822.1] from China. Anopheles sp from Ratahan based on the COI gene has a similarity level of 
6096 with Anopheles barbirostris [KM610029.1] from China. Anopheles sp from Pineleng has a 77% 
similarity with Anopheles aquasalis [AF417697.1] from Brazil. The COI gene sequences of Anopheles sp 
from Southeast Minahasa (Tombatu and Ratahan), and Anopheles sp from Minahasa (Pineleng) had a 
nitrogen base size difference of more than 6%. Thus, the variation of the Anopheles sp COI gene is 
relatively high compared to similar sequences that have been recorded on the NCBI gene bank site. 


Keywords— Reconstruction of Phylogeny, Anopheles sp. Cytochrome Oxidase Sub Unit 1 (CO1) gene, 
Minahasa 


I. INTRODUCTION endemic areas in Indonesia. Some malaria endemic areas 
Malaria is still a major health consideration, especially in in North Sulawesi are Minahasa, Southeast Minahasa and 


tropical countries. Malaria is the world's most dangerous North Minahasa. Malaria is an infectious disease caused by 


parasitic infection, causing more than a million death and a protozoan parasite of the genus Plasmodium, which is 


500 million cases annualy (Penet et al.,2007; Ravichandran transmitted through the bite of the Anopheles mosquito. 


et.al. 2007). Malaria may decrease the productivity of Malaria can only be transmitted through the bite of the 


individuals, families and the whole due morbidity and Anopheles Sp. Thus, the high case of malaria infection in 


mortality (Ravichandran et.al. 2007, Namdeo, et.al., 2006). an area 1s linear with the high population of Anopheles sp. 


Malaria remains a leading cause of morbidity and mortality In Indonesia, vector confirmation has been carried out 


worldwide with an estimated 500 million cases and 2.5 from 1919 to 2009, and during that period 25 species were 


million deaths annually (Stauffer & Kamat 2003). Malaria found to be positive for the malaria parasite. As a tropical 
is a reemerging disease, which is a disease that is re- rain forest area, the Minahasa area is separated by forests 
and mountains. Based on a survey from the Ministry of 
Health of the Republic of Indonesia in 2009, in North 


Sulawesi, three main species of Plasmodium vector were 


infected en masse (Arsin, 2012). 
Indonesia is a country with the highest malaria cases in the 
world. North Sulawesi is known as one of the malaria 
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found in humans, namely Anopheles subpictus, Anopheles 
vagus and Anopheles annullaris (Ministry of Health RI, 
2009). Identification is based on morphological 
characteristics. 

However, identification was carried out based on 
morphological characteristics. Identification of mosquitoes 
by morphological analysis method has many limitations. 
The observed specimens often have undergone 
morphological changes due to immersion with alcohol or 
formalin from the sampling location. This greatly affects 
the stage of species identification in the laboratory. 

The identification method to the species level that has high 
accuracy is by molecular identification. Molecular 
identification using the cytochrome oxidase sub unit 1 
(CO1) gene has been widely carried out. CO1 is one of the 
genes in mitochondrial DNA that has been designated as a 
molecular barcode. The cytochrome C oxidase sub unit I 
(COLI) gene has special characteristics that are suitable as a 
tool in evolutionary studies, namely (1) as a final catalyzer 
in the respiratory chain in mitochondria, so COI is widely 
studied at the biochemical level, and shows that the 
structure and size of the COI gene conserved in all aerobic 
organisms (Rivera et al. 2009). (2) The amino acid 
sequence correlates with the function of each part of the 
COI, thus showing the characteristics of the species that 
possess it (Rivera et al. 2009 Roe & Sperling 2007). (3) A 
sequence of 658 basepair (bp) at the 5' end was proposed 
as an animal barcode (Hebert et al. 2003 a,b). These 
barcodes have been successfully demonstrated to be able to 
differentiate between species in Lepidoptera (Hebert et al. 
2003a; Hajibabaei et al. 2005), beetles (Funk et al. 1995), 
several insect pests (Toda & Murai 2007) moth Hamona 
mermerodes (Hulrc et al. 2007). al. 2007), mosquitoes 
(Cywinska et al. 2006). 

Identification of insects from North Sulawesi using the 
COI gene has been successfully carried out on Aedes sp 
(Kaunang etal. 2013; Timah et.al. 2016), Apis dorsata 
Binghami (Mokosuli et.al. 2013), subterranean termites 
(Ngangi et.al. 2014), Demselfly (Rantung et.al. 2015) and 
Marine Gerridae (Warouw et.al. 2015). Molecular 
identification using the COI gene as the basis for 
reconstructing the phylogeny of Aedes sp. in North 
Sulawesi as a malaria endemic area. Reconstruction of 
phylogeny will break the distribution of Anopheles sp. in 
Minahasa. The results of the phylogeny reconstruction will 
be very useful for the prevention of malaria vector 
mosquitoes. 


II. MATERIALS AND MEHODS 
Sample 
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Adult mosquito collected used modified method Cheng et. 
al. (2010). Collection on the fields area randomly. Insects 
that have been collected will insert in a bottle sample that 
has been labeled with place and time of data sampling. The 
bottle was filled with 70% alcohol for identification and 
preservation. 


DNA Extraction, PCR Amplification and Sequencing 


Total genomic DNA was extracted larva and adults 
mosquito using Qiagen DNA Blood and Tissue, according 
to the manufacturer’s protocol. PCR was performed in a 
total volume of 25 uL containing 1 x reaction buffer, 3 
mM MgCl2, 0.24 mM dNTPs, 1.4 uM of each primer 
LCO1490 : 5-GGTCAACAAATCATAAAGATATTGG- 
3’ and HCO2198 5'- 
TAAACTTCAGGGTGACCAAAAAATCA-3'  (Folmer 
et. al, 1994), 1U Go Taq Flexi DNA polymerase 
(Promega Corp.) and 2.5 uL of DNA (a 100 time dilution 
of the original DNA). The PCR program was as follows: 
94 °C for 5 min, followed by 40 cycles of 94 °C for 1 min, 
48 °C for 1 min and 72 °C for 1 min and a final extension 
at 72 °C for 5 min. PCR products were purified using 
Wizard SV Gel and PCR Clean-Up System (Promega 
Corp). Purified PCR products were analyzed by 
electrophoresis in 1% agarose gel. The molecular size of 
the amplified products was estimated using 1 kbp DNA 
ladder (Biometra). PCR products were sequenced using 
ABI PRISM Dye Terminator Cycle Sequencing Ready 
Reaction System, version 1.1. (Applied Biosystems) in 
FIRST BASE Singapura 


Sequences Analyses and Phylogeny trees reconstructed 


Obtained sequences were aligned using MEGA 6.0 and 
Geneous 6.0 software. Sequences were subjected to Basic 
Local Alignment Search Tool (BLAST) in order to 
perform sequence similarity searches 
(www.ncbi.nih.gov.com). Nucleotide frequencies were 
calculated using MEGA 6.0 software (Tamura et. al. 
2013). The genetic distances (number of nucleotide 
substitutions per site) among sequences were calculated 
using the Maximum Composite Likelihood model in 
Geneous 6.0 software. Phylogenetic trees were 
reconstructed using two different reconstruction methods: 
(1) neighbor joining (NJ) and (2) maximum parsimony 
(MP). The NJ tree was reconstructed using the Maximum 
Composite Likelihood method. Phylogenetic analyses 
were conducted in MEGA 6.0 software. Bootstrap support 
values were obtained by 1,000 replications using both 


methods (Tamura et. al. 2013). 


III. RESULTS AND DISCUSSION 
Total DNA Extraction and Purification 
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Tissue in adult insects is found in the exoskeleton. Thus, to 
obtain good DNA purity and concentration, the selection 
of the right organ will determine the success of DNA 
extraction. In this study, a series of trials were carried out 
using mosquito organs to obtain total dsDNA with 
standard purity and concentration. The total DNA in this 
study was nuclear DNA and mitochondrial DNA isolated 
from mosquito organ cells which were extracted using the 
DNA blood and tissue kit. The organs used are the head, 
thorax, legs and abdomen. From the extraction carried out, 
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it was found that the use of the thoracic organ resulted in 
the best concentration and purity (table 1). Proteinase K is 
a key enzyme in DNA extraction. Proteinase-K functions 
to damage membrane proteins and other proteins in the 
cytoplasm and nucleoplasm as well as other cell 
compartments to isolate DNA. The results showed that 
proteinase K immersion time affected the total DNA purity 
but had no effect on the total DNA concentration obtained 
(Table 1). 


Table I. Purity and Concentration of total Mosquito DNA 


No Ortum Sample weight Modified Immersion Purity Concentration 

(mg) Proteinase-K (hours) (A260/A280) ug/ml 

1 Head 35 12 1,10 25,80 
35 24 1,23 25,84 

2 Thorax 35 12 1,72 37,50 
35 24 1,87 36,24 

3 Limbs 35 12 1,35 32,10 
35 24 1,45 33,40 

4 Abdomen 35 12 1,42 34,50 
35 24 1,52 36,20 


Amplification and Visualization of Gene CO1 


The results of the CO1 gene amplification were shown by 
electrophoresis electrogram. Based on the formed band, it 
shows a high concentration of amplicon in both samples 
ARI, KLI and SG1 (Figure 1) 





Fig.l Visualization of the PCR product of Anopheles sp 
COI] gene amplicons through 0.596 agarose gel 
electrophoresis. Anopheles sp from Pineleng/PSA (no. 9), 
Anopheles sp from Ratahan/RTA (no. 10) and Anopheles 
sp from Tombatu/TLA (no. 11). 
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The results of partial sequencing of Gen CO1 in the form 
of an ABI file, interpreted using Geneous 6.0 software. 
The sequence lengths of RTA (Ratahan), TLA (Tombatu) 
and PSA (Pineleng) were 870 bp, 774 bp and 862 bp, 
respectively. Based on the sequencing chromatogram, it 
showed that the CO1 gene sequence formed was good 
(Figure 2). 
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Fig.2c. The nitrogen base sequence of the PSA COI gene was read with the Geneous Program 6.0 
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Table 2. Percentage similarity of TLA CO1 gene sequences compared with the top ten Sequences recorded in the NCBI gene 
bank (https://blast.ncbi.nlm.nih. gov/Blast.cgi) 
N 
4 


f 


1 Anopheles maculatus voucher 9 KT382822.1 
AMAC20150811V4 mitochondrion, 
complete genome 
Anopheles albitarsis cytochrome c oxidase AF417696.1 
subunit I (COI) gene, partial cds; 
mitochondrial gene for mitochondrial 


product 
DQ076230.1 


Anopheles deaneorum mitochondrion, HQ335347.1 
complete genome 
5 Anopheles marajoara cytochrome c oxidase AF417699.1 
subunit I (COI) gene, partial cds; 
mitochondrial gene for mitochondrial 
product 


Anopheles deaneorum isolate D6 
cytochrome oxidase subunit I (COI) gene, 
partial cds; mitochondrial 


Anopheles albitarsis isolate A3 cytochrome DQ076206.1 
oxidase subunit I (COD) gene, partial cds; 
mitochondrial 


Anopheles aquasalis isolate GUA109012 


OO 
OO 


KC354821.1 


cytochrome oxidase subunit I (COI) gene, 


partial cds; mitochondrial 


Anopheles oswaldoi haplotype H11 
cytochrome oxidase subunit I gene, partial 
cds; mitochondrial 


OO 
OO 


DQ784837.1 


OO 
OO 


OO oo oo OO 
OO oo oo oo 


DQ784836.1 


Anopheles oswaldoi haplotype H10 


o 
2 
3 
7 
cytochrome oxidase subunit I gene, partial 
cds; mitochondrial 
10 Anopheles marajoara isolate C8 DQ076223.1 
cytochrome oxidase subunit I (COI) gene, 


partial cds; mitochondrial 
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Table 3. Percentage similarity of COI RTA gene sequences compared with the top ten Sequences recorded in the NCBI gene 
bank (https://blast.ncbi.nlm.nih. gov/Blast.cgi) 


No 
1 Anopheles aquasalis isolate aqua28 4e137 AF548901.1 
cytochrome oxidase subunit I (COI) gene, 
partial cds; mitochondrial gene for 
mitochondrial product 
2 Anopheles aquasalis isolate aqua10 4e137 
cytochrome oxidase subunit I (COI) gene, 
partial cds; mitochondrial gene for 
mitochondrial product 
3 
5 
7 
10 


E AF548894.1 
Anopheles barbirostris subgroup clade III 5e136 83 KM610037.1 
isolate SMMULZ3 
cytochrome oxidase 
subunit I (COI) gene, partial cds; 
mitochondrial 
Anopheles barbirostris subgroup clade III 5e136 83 KM610022.1 
isolate SMMUPR3 
cytochrome oxidase 
subunit I (COI) gene, partial cds; 
mitochondrial 
Anopheles oswaldoi haplotype H09 5e136 83 DQ784835.1 
cytochrome oxidase subunit I gene, partial 
cds; mitochondrial 
E 
- 
A i 


4 


DQ784834.1 


Anopheles oswaldoi haplotype H08 5e136 
cytochrome oxidase subunit I gene, partial 
cds; mitochondrial 


Anopheles barbirostris subgroup clade II 4e142 
isolate SMMUPR10 

cytochrome oxidase 

subunit I (COI) gene, partial cds; 

mitochondrial 


Anopheles barbirostris subgroup clade III 6e140 
isolate th1.10 cytochrome oxidase subunit I 
(COD) gene, partial cds; mitochondrial 


Anopheles barbirostris subgroup clade III 1e136 


KM610029.1 


EU797223.1 





EU797224.1 


isolate kh3 cytochrome oxidase subunit I 


(COI) gene, partial cds; mitochondrial 





Anopheles barbirostris subgroup clade III 1e136 EU797218.1 
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isolate th1.9 cytochrome oxidase subunit I 


(COI) gene, partial cds; mitochondrial 
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Table 4. Percentage similarity of PSA COI gene sequences compared with the top ten sequences recorded in the NCBI gene 
bank (https://blast.ncbi.nlm.nih. gov/Blast.cgi) 


Anopheles aquasalis isolate aqualO 
cytochrome oxidase subunit I (COI) gene, 
partial cds; mitochondrial gene for 
mitochondrial product 

Anopheles aquasalis isolate aqua21 
cytochrome oxidase subunit I (COI) gene, 
partial cds; mitochondrial gene for 
mitochondrial product 

Anopheles aquasalis cytochrome c oxidase 
subunit I (COI) gene, partial cds; 
mitochondrial gene for mitochondrial 
product 

Anopheles oswaldoi haplotype H05 
cytochrome oxidase subunit I gene, partial 


cds; mitochondrial 


5 Anopheles oswaldoi haplotype H01 
cytochrome oxidase subunit I gene, partial 
cds; mitochondrial 


Anopheles oswaldoi haplotype H04 


cytochrome oxidase subunit I gene, partial 
cds; mitochondrial 

Anopheles oswaldoi haplotype H02 
cytochrome oxidase subunit I gene, partial 
cds; mitochondrial 

Anopheles oswaldoi haplotype H03 
cytochrome oxidase subunit I gene, partial 
cds; mitochondrial 

Anopheles marajoara isolate C1 
cytochrome oxidase subunit I (COI) gene, 
partial cds; mitochondrial 

Anopheles punctipennis cytochrome c 
oxidase subunit I (COI) gene, partial cds; 
mitochondrial gene for mitochondrial 


product 


ISSN: 2456-1878 (Int. J. Environ. Agric. Biotech.) 
https://dx.doi.org/10.22161/ijeab.65.10 





2e100 78 


7e100 


9e104 


AF548894.1 


oo 


a a - Q T - Q ~ - 
ON ON 3 3 XQ 3 J - 
Ui 


AF548900. 1 


AF417697.1 


DQ784831.1 


DQ784827.1 


DQ784830.1 


DQ784828.1 


DQ784829.1 


DQ076216.1 


AF417720.1 


Marthy Lingkan Stella Taulu et al. 


The results of BLAST NCBI analysis of the COI gene 
sequences of Anopheles sp from Tombatu showed a 92% 
similarity with Anopheles maculatus [KT382822.1] from 
China. Anopheles sp from Ratahan based on the COI 
gene has a similarity level of 80% with Anopheles 
barbirostris [KM610029.1] from China. Anopheles sp 
from Pineleng has a 7796 similarity with Anopheles 
aquasalis [AF417697.1] from Brazil. (Table 4, Table 5 
and Table 6). Alignment results of Anopheles sp from 
Tombatu showed 6 different nitrogen base sites with 
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sequences similar to Anopheles maculatus [KT382822.1]. 
Meanwhile, Anopheles sp from Ratahan showed 39 
different nitrogen base sites with similar sequences to 
Anopheles barbirostris [KM610029.1]. Anopheles sp 
from Pineleng showed 38 different nitrogen base sites 
with sequences similar to Anopheles  aquasalis 
[AF417697.1] (Table 4, Table 5 and Table 6). The 
position of the difference in nitrogen bases indicates a 
mutation that occurs in Anopheles sp from Minahasa. 


Table 5. Alignment of Anopheles from Pineleng with Sequences Similar to NCBI Anopheles aquasalis. 
(https://blast.ncbi.nlm.nih.gov/Blast.cgi) 


Anopheles aquasalis cytochrome c oxidase subunit | (COI) gene, partial cds; mitochondrial gene for mitochondrial product 


Sequence ID: AF417697.1 Length: 899 Number of Matches: 1 


Range 1: 132 to 884 


Score Expect Identities Gaps Strand Frame 
448 bits(242) 18-1210) 595/TT AT 796) INT 3(496) Plus/Minus 
Features: 


seset eoe LE CT UHL es 
ee 35 amem ENS 


Sbjct 824 


765 


rr c TITE TINT n 


Sbjct 764 


err nu ME TT TEE [RTT beni pape os 


Sbjct 704 GTCGAGGTATT 
Query 275 N-TCHN-CAAATATTATATCAAATTGT 


705 


645 
TTTAAMNATGTTGG-TTTÀ 331 


icc eaa lad chen MHA henL MH chalets see 
rr es OT HIT IT TTT TP Tr Ty Te 


Sbjct 585 TCCTGTTAAAAGAGG 


CTATAATAGCAAATAC- 528 


wey se SETTTTITHeT APATE TIPP PTHTP A HT TT nn 


Sbjct 527 CCTATTGATAATACATAATGGAAATGA-GCTACTAC-ATAATATGTGTCATGTAGT-ACA 471 


e ET TTT TU Pir Los 


Sbjct 478  ATGTCAATTGAAGAATT, 
Query 58 


j 


CCCAGTTAATCCACCTACAGT -À 413 


VITRO TETTE METETE een Te 2 


Sbjct 412 IVA RUUTON -AAAGTATAGCTGGGCTAT-A-TGTTAATTGTGTTCCATG | 358 


er HEAT TENET n 
ce or ETE ANE 
ce or EON 7 
e cm TRIER c 
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Table 6. Alignment of Anopheles from Tombatu with Similar Sequences on NCBI Anopheles aquasalis. 
(https://blast.ncbi.nlm.nih.gov/Blast.cgi) 
Anopheles maculatus voucher AMAC20150811V4 mitochondrion, complete genome 
Sequence ID: KT382822.1 Length: 14850 Number of Matches: 1 
Range 1: 2177 to 2818 


Score Expect Identities Gaps Strand Frame 
1044 bits(565) 0.0() 686/745(9296) 6/745(0% ) Plus/Minus 
Features: 


me aso MULTUM HEEL A RH rA LAL zoe 


ia sese err rir RARE 28 
Query 356 AATTEENA TTAN 20 


E A A ATA As ase 
sje aer AULAE HH. asa 
sje 2 DAMALS ase 


Query 398 CATAATGAAAAT(G 449 


sosce asso dads LULA aso 
cory TI TTC 2 
seset 2 UBL A AAA HR 220 
pera PH RT 2520 


"| 


al 


ve zuo. HL eA Led cjn dM 226 


e c EHTEL TTT. 
er e Cesare age 7 


bjct 2288 n AAA ATE- AGAAAT 
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Table 7.Alignment of Anopheles of Ratahan with Similar Sequences in NCBI Anopheles aquasalis. 
(https://blast.ncbi.nlm.nih.gov/Blast.cgi) 


Anopheles barbirostris subgroup clade lll isolate th1.10 cytochrome oxidase subunit | (COI gene, partial eds; mitochondrial 
Sequence ID: EUT97223.1 Length: 756 Number of Matches: 1 


Range 1: 62 to 756 

















aoe as tir SEDENT D o 
e OE l| 2. 
esc cor MOLLE CULO TTT mmm " 
eser eo dL EAT TERR ecl TE aea 
scr ass Dallas EE MAL Ae EMI AD TU aes 
cr aon Asc ELA En ue ETT aas 


This genetic variation is supported by the results of 
morphometric analysis which indeed show differences in 
several morphometric characters, including the shape of 
the strip on the pronotum and the color of the antennae 
tip on Rhynchoporus sp. cream-colored sago palm and 
Rhynchoporus sp. in black sugar palm (Korua et.al. 
2015). Polymophism can occur in a population if more 
than one morphological variation is found at the same 
location and time (Ford, 1965, Abad et. al. 2014). If 
random mating occurs and each individual has the 
potential to mate, then morphological changes can take 
place in a population (Abad et. al. 2014). 


Phylogenetic Analysis and Construction 


The substitution matrix between Anopheles sp from 
Tombatu, Ratahan and Pineleng was compared with 22 
BLAST sequences at the NCBI site built using the 
Maximum Likelihood Model on the MEGA 6.0 program. 
The form of transitional substitution is indicated by the 
numbers in bold in table 6. While the transversional 
substitution is written in italics in table 6. Nucleotide 
frequency A = 33.70%, T/U = 36.75%, C = 14.11% and 
G = 15.44 %. The maximum value of the Log Likelihood 
from the calculation results is 9678, 156. 
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Table 6. Estimation of the Maximum Likelihood Model 
Substitution Matrix in the MEGA 6.0 . Program 


A T/U C G 
A - 8.24 3.16 3.71 
T/U 7.55 12.03 3.46 
C /.923 31.33 - 3.46 
G 8.09 8.24 3.16 


Phylogenetic construction 


The phylogeny construction was carried out using two 
models, namely Neighbor Joining and UPGMA. These 
two models are used because they have similarities, 
namely the evolutionary approach and to see the position 
of the Anopheles sp species from Minahasa. The 
Neighbor Joining model phylogeny construction was 
built with 22 sequences similar to the NCBI BLAST 
results. Three monophyletic clades were formed, where 
Anopheles sp from Tombatu and Ratahan were in the 
same clade while Anopheles from Pineleng was in its 
own clade. The phylogenetic tree construction using the 
UPGMA model also placed Anopheles sp from Pineleng 
in its own node, while Anopheles sp from Tombatu and 
Ratahan formed the same node but stil in one 
monophyletic clade. In the UPGMA model, only 2 
monophyletic clades were formed. 
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DQ076206.1| Anopheles albitarsis 
DQ076223.1| Anopheles marajoara 
DQ076230.1| Anopheles deaneorum 
DQ784835.1| Anopheles oswaldoi 
DQ784837.1| Anopheles oswaldoi 
DQ784836.1| Anopheles oswaldoi 
KT382822.1 Anopheles maculatus 
HQ335347.1| Anopheles deaneorum 
KT382822.1| Anopheles maculatus 
Ratahan 
Tombatu 
AF417696.1| Anopheles albitarsis 
AF417697.1| Anopheles aquasalis 
AF417699.1| Anopheles marajoara 
AJ271384.1| Anopheles dirus 
KC354821.1| Anopheles aquasalis 
EU797223.1| Anopheles barbirostris 
EU797224.1| Anopheles barbirostris 
EU797218.1| Anopheles barbirostris 
KM610029.1| Anopheles barbirostris 
KM610037.1| Anopheles barbirostris 
KM610022.1| Anopheles barbirostris 
Pineleng 
AF548901.1| Anopheles aquasalis 
AF548894. 1| Anopheles aquasalis 


Fig.3: Anopheles sp phylogeny tree from Tombatu, Ratahan and Pineleng compared to 22 BLAST sequences at the NCBI 
site, built using the Neighbor Joining Model, bootstrap 1000 x. 


DQ076230.1| Anopheles deaneorum 
DQ076223.1| Anopheles marajoara 
DQ076206.1| Anopheles albitarsis 
DQ784835.1| Anopheles oswaldoi 
DQ784837.1| Anopheles oswaldoi 
DQ784836.1| Anopheles oswaldoi 
HQ335347.1| Anopheles deaneorum 
KT382822.1 Anopheles maculatus 
KT382822.1| Anopheles maculatus 
Pineleng 

AF548901.1| Anopheles aquasalis 
AF548894.1| Anopheles aquasalis 
Ratahan 

Tombatu 

EU797224.1| Anopheles barbirostris 
EU797218.1| Anopheles barbirostris 
EU797223.1| Anopheles barbirostris 
KC354821.1| Anopheles aquasalis 
AF417699.1| Anopheles marajoara 
AJ271384.1| Anopheles dirus 
AF417696.1| Anopheles albitarsis 
AF417697.1| Anopheles aquasalis 
KM610022.1| Anopheles barbirostris 
KM610029.1| Anopheles barbirostris 
KM610037.1| Anopheles barbirostris 


Fig.4: Anopheles sp phylogenetic tree from Tombatu, Ratahan and Pineleng compared to 22 BLAST sequences at the 
NCBI site, constructed using the UPGMA Model, bootstrap 1000 x. 
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IV. CONCLUSION 


1. Based on the COI gene, the mitochondrial DNA of 
Anopheles sp from Tombatu is 92% similar to Anopheles 
maculatus [KT382822.1] from China. Anopheles sp from 
Ratahan based on the COI gene has a similarity level of 
80% with Anopheles barbirostris [KM610029.1] from 
China. Anopheles sp from Pineleng has a 77% similarity 
with Anopheles aquasalis [AF417697.1] from Brazil. 


2. The COI gene sequences of Anopheles sp from 
Southeast Minahasa (Tombatu and Ratahan), and 
Anopheles sp from Minahasa (Pineleng) have nitrogen 
base size differences of more than 6%. Thus, the variation 
of the Anopheles sp COI gene is relatively high 
compared to similar sequences that have been recorded on 
the NCBI gene bank site. 
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